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Abstract 

It is demonstrated how perturbative 0(a 2 ) relativistic corrections to the non-relativistic stable 
heavy quark-antiquark production cross section in e + e~ annihilation based on a Coulombic QCD 
potential can be systematically calculated using the concept of effective field theories. The 0{a 2 s ) 
corrections from the relativistic energy-momentum relation, the relativistic phase space corrections 
and the 0{a 2 s ) corrections involving the group theoretical factors C\ and CfT are determined ex- 
plicitly. For the case of ti production the sum of these corrections amounts to 3% - 7% over the 
whole threshold regime and is insensitive to variations in the top width. Perturbative 0(a 2 s ) cor- 
rections to the leptonic decay width of heavy quark-antiquark 3 S\ vector resonances are extracted. 




In the kinematic regime where the c.m. energy is much larger than the quark masses the total hadronic 
cross section in e + e~ collisions Rhad = o"(e + e~ — ► Hadrons)/<7 p t, where a p t is the point cross section, 
belongs to the theoretically best known and understood quantities in electroweak physics. Inspired 
by the asthonishing precision of the LEP experiments Rhad has been calculated up to 0(a^) accuracy 
for LEP energies based on sophisticated multi-loop techniques [||. In the kinematic regime near 
quark-antiquark thresholds, however, which is characterized by the relation 

V s + te 

where Mq is the quark mass^j and yfs the c.m. energy, the theoretical understanding of Rhad is 
much poorer. Although perturbative methods to describe the hadronic cross section in the threshold 
regime in general seem to be more than questionable due to non-perturbative effects and due to the 
large size of the strong coupling, it is well accepted that a perturbation theory based description 
of the hadronic cross section near a quark-antiquark threshold should be possible if the quark mass 
is sufficiently larger than the hadronic scale A QCD [0]. The prototype application of perturbative 
methods in the threshold regime is therefore the case of ti production^ since non-perturbative effects 
are suppressed due to the large top quark mass and width || U and since the strong coupling is 
sufficiently small, a s (CFMta s ) ~ 0.15. The most common approach found in literature is based on 
numerical solutions of a non-relativistic Schrodinger equation with a phenomenological QCD potential 
where the short-distance part originates from loop calculations and the long-distance part is obtained 
from fits to charmonium and bottomonium spectra 0. Although this approach leads to applicable 
results, the use of a phenomenological potential makes the systematic implementation of dynamic as 
well as kinematic relativistic corrections difficult, if not impossible. As a consequence all potential 
model predictions for the cross section in the threshold region contain intrinsic uncertainties of relative 
order \(3 2 \ a 2 . The actual size of these 0(a 2 ) corrections^] and their dependence on the c.m. energy 
is unknown because a consistent and systematic determination of these corrections has never been 
achieved. It is the purpose of this letter to demonstrate how relativistic perturbative 0(a 2 ) QCD 
corrections to the non-relativistic stable heavy quark-antiquark cross section predicted by the non- 
relativistic Schrodinger equation with a Coulombic QCD potential V^cd = —Cpa s /r and a fixed a s 
can be calculated. We want to emphasize that this work is not meant to present all calculational 
details, but to show the main steps and results of our calculations. A more detailed and explicit work 
will be published later. We also would like to mention that at no point in this work non-perturbative 
effects and electroweak corrections are taken into account. 

The formula for the single photon mediated (i.e. vector-current induced) non-relativistic cross 
section valid also for complex energies reads S 



rQQ, 



a{e + e~ QQ) 

(Tpt 



7+ln(-i/3)+*(l- 



1 Throughout this paper the quark mass Mq is understood as the pole mass and /3 is called "velocity" . 

2 An application of perturbative methods might also be possible for the bb production cross section for low radial 
excitations of the T family because in that case non-perturbative effects seem to be well under control JE], ^ and the 
strong coupling is still small enough that a perturbative calculation might be justified. 

3 Relativistic corrections to the non-relativistic cross section actually result in a series in powers of a s and f3. Because 
we consider the kinematic regime \/3\ < a s , we count powers of fi as powers of a s . For simplicity we will call the C(a 2 ), 
0{a s j3) and 0(f3 2 ) relativistic corrections briefly "C(a 2 ) corrections". 
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where is the digamma function, ^(z) = 4-lnT(z), <r p t = 4-7ra 2 /(3s) and A^ c = 3 is the number 
of colors, eg is the electric heavy quark charge and a denotes the fine structure constant. In the 
non-relativistic limit and for unstable quarks, where the relation between the energy E = y/s — 2Mq 
relative to the threshold point and the velocity can be approximated as j3 = [{E + iTq)/Mq\ 1 / 2 , Tq 
being the width of the heavy quark, formula @ coincides with an expression given in [||, Q. For stable 
quarks and above threshold expression (|2|) leads to the famous Sommerfeld expression 

r qq,nr _ 3 2 C F a s -K 

whereas for stable quarks and below threshold eq. @ develops narrow resonances at the well known 
Coulomb energy levels [0, M. So far only 0(a s ) QCD corrections to the non-relativistic cross section 
from the running of the strong coupling |7], |l^, |5| and from short distances have been successfully 
calculated. The knowledge of the 0(a 2 ) corrections to expressions @ and (||) is important in order 
to understand the structure and importance of relativistic corrections and might even lead to hints on 
the form of non-perturbative effects. Because it is quite unclear how the 0(a 2 ) relativistic corrections 
to the non-relativistic cross section obtained from a Coulombic potential have to be implemented 
consistently into a potential model approach they will serve as an order of magnitude estimate for the 
uncertainties inherent in potential model predictions. 

In this letter we explicitly determine the 0(a 2 ) corrections to the total single photon medi- 
ated heavy quark pair production cross section from the relativistic energy-momentum relation, the 
relativistic phase space corrections and those 0(a 2 ) corrections which are multiplied by the SU(3) 
group theoretical factors C F and CfT (Cf = 4/3, T = 1/2), where for the latter contribution only 
the effects from the heavy quark itself are taken into account. We would like to stress that these 
corrections represent a gauge invariant subset of all 0(a 2 ) QCD corrections and that no model- like 
assumptions are imposed for our calculation. The 0(a 2 ) QCD corrections involving the SU(3) group 
theoretical factors CaCf and CfTui, where Ca = 3 and ni is the number of massless quark flavors, 
including the 0(a 2 ) effects from the running of the strong coupling are not treated in this work, but 
can be determined along the lines presented here. For the case of a non-negligible width of the heavy 
quark (as in ti production) further types of corrections have to be taken into account. Apart from 
the perturbative (multi-loop) corrections to the width of the free heavy quark also corrections from 
the off-shellness of the decaying heavy quark, from the interaction among the decay products and the 
other heavy quark (if it is not decayed yet) and from time dilatation effects need to be considered. 
Although the size and the interplay of all these effects have been examined at various places in liter- 
ature (see e.g. [12, O, 14]) they are still not completely understood yet as far as 0(a 2 ) corrections to 



the cross section are concerned. However, it is agreed in [12, [b|, 14 that the bulk of these effects can 
be accounted for by using a momentum dependent width instead of a constant one. This difficulty is 
completely ignored in this work. Therefore, all formulae presented here are strictly speaking only valid 
for stable quarks. However, during our calculations we never assume that the squared velocity f3 2 is a 
real number. Therefore, our results could be easily implemented in a more complete approach which 
treats all finite width effects properly. For now we will use the naive replacement E — > E + i Tq in the 
spirit of j|, H|, where Tq represents an appropriately chosen constant heavy quark width, which is not 
necessarily equal to the decay width of a free heavy quark. We think that this procedure is justified 
in order to demonstrate the size of the 0(a 2 ) corrections calculated in this work in the presence of a 
large quark width. We finally would like to emphasize that relativistic corrections from the exchange 
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of non-instantaneous gluons (responsible for Lamb shift type corrections and real radiation effects for 
s > 4Mq) are of order |i5| , |i~6[ ] and beyond the level of accuracy intended in this work. 

Let us start by reminding the reader that by means of the optical theorem the non-relativistic 
cross section ^QQ> NR ; e q. (H) 5 can be written as [||, [|] 



Im 



N r e, 



2 2 -G c (0,0) 



(4) 



where G c is the Coulomb Green function satisfying the equation of motion 



1 - 2 a p 2 
Mq~ g ~W\~MQ 



G c (x,y) = 6®(x-y). 



a = Cfol s 



(5) 



p is called the "external three momentum" for the rest of this work. It should be noted that p 2 /Mq is 
equal to the energy E = ^fs — 2Mq only up to higher orders in E. In quantum mechanics textbooks 
usually the convention is employed where these higher order terms are set to zero. As explained later, 
we set p 2 equal to the squared (relativistic) three momentum of heavy quarks in the cm. frame. 
An explicit analytic expression for G c in coordinate space representation has to our knowledge been 
calculated the first time in |l7j. Our strategy is to calculate relativistic corrections to the Coulomb 
Green function G c by using (textbook quantum mechanics) time independent perturbation theory 
(TIPT). This approach is justified because only the exchange of instantaneous gluons needs to be taken 
into account. Using the concept of effective field theories the arising UV divergences are removed by 
matching to recent two-loop results for the cross section in the threshold region in the framework of 
QED [15]. The cross section is then obtained via the optical theorem, eq. (||). Conceptually we follow 
the lines presented in ||. 

The relativistic corrections to the Coulomb Green function arise from three different sources: 
(a) the relativistic energy-momentum relation, (b) 1/Mq corrections to the Coulomb potential Vq CD 
and (c) 1 /Mq corrections from the electromagnetic current which produces and annihilates the quark- 
antiquark pair. The corrections from the relativistic energy-momentum relation can be easily deter- 
mined by taking into account that G c can be written in the form M 



G c (x,0) 



d 6 p eti 



(2vr) 3 



vpox 



Mr 



e n /,/V " 



An a 



Mr 



n=i J ( 27r ) 3 (Pn-1 ~ Pn) 2 Pn ~ P 2 + ^ 



(6) 



The relativistic expression for the particle-antiparticle propagation J dp^ / (2ir) 3 Mq / (p 2 — p 2 + ie) 
reads 

~i j^A S {Pn+ (.JM 2 + p 2 , ) ) S(pn ~ (.JM 2 + p 2 , ) ) , Pn (?) 

where S is the common Dirac propagator. Because the Coulomb interaction is instantaneous the 
Pn integration can be carried out. An expansion in 1/Mq then yields that the corrections from the 
relativistic energy-momentum relation can be implemented into expression (H) by the replacement 



Mq 



Mq 



Pn — P 2 + i e 



Pn — P 2 + it 



1 + 



Pn+P 2 

AM 2 



(8) 



for each particle-antiparticle propagator. As mentioned earlier, the form of the correction factor on the 
r.h.s. of eq. (H) differs from the usual kinetic energy correction used in quantum mechanics textbooks 
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because the relation between the external momentum p and the cm. energy reads 



M 2 



1-P 



2 ' 



(9) 



For stable quarks definition @ leads to the relation p 2 /Mq = E + E 2 /(4Mq) + 0(E 3 ) between the 
external momentum p and the energy E. Any other definition of p 2 would lead to a different form 
of the r.h.s. of eq. (|8|). The final result for the cross section, of course, is independent of this choice. 
We have chosen definition (^) to facilitate our calculations. The 1/Mq corrections to the interaction 
potential are well known and read 



V(Q) 



Air a ir a Air a 



Q 2 

47TO 



Mr 



P 2 



Mr 



S\ S2 



(QS 1 )(QS 2 ) 



(OP 
Q 4 



\2 1 



6ira 

~mJ 



a 2 



(Si + s 2 [ 



Q x P 



(10) 



in momentum space representation, where Q is the (three) momentum flowing through the gluons and 
the Su2 represent the quark/antiquark spin operators. In eq. (|l0|) the Coulomb interaction is also 
displayed. The 1/Mq corrections from the electromagnetic current which produces and annihilates the 

-*2 — * — * 

heavy quark-antiquark pair lead to the insertion of the factor {1 — ^r[| + >Si<%]} into expression (g). 

Taking into account that the QQ pair is produced in a (J PC = 1 ) 3 Si state^the relativistic 
corrections to the Coulomb Green function can be rewritten in terms of corrections induced by an 
effective interaction potential which in coordinate space representation takes the simple form 



V 3 si(x) 



a 

\x\ 



1 + 



3p 2 1 

2M} 



11 Tra_ (3) 
+ 3 M 2 ° [X) 



4 Mr 



(11) 



where the Coulomb potential is also displayed. In addition, there remains a correction to the Coulomb 
Green function, which cannot be expressed in terms of an interaction potential. This correction takes 
the form 



lim 

|x|-»0 



1 



Mr 



7 -v 2 +p 2 

o 



G c {x,Q) 



(12) 



and essentially represents relativistic phase space corrections. 

It is now straightforward to determine all 0{a 2 ) corrections to the Coulomb Green function. 
The corrections from the first term on the r.h.s. of (|Tl"D, called kinetic energy corrections later in this 
paper, can be trivially implemented by the replacement a — * a(l + 3p 2 /2Mq) in G c . The corrections 
from the second and third term, later called dynamical corrections, can be calculated via coordinate 
space TIPT. The phase space corrections of eq. (|TJ) can be evaluated by employing the equation of 
motion d). The arising UV (short-distance) divergences can be regularized by considering corrections 
to G c (x, 0) and taking the limit \x\ — > afterwards^. Taking into account that the relation between 

4 The production of a 3 Di state is proportional to the modulus squared of the second derivative of the heavy quark- 
antiquark wave function at the origin and therefore suppressed by |/3| 4 w a 4 . This is beyond the intended accuracy. 

5 We emphasize that this method is not claimed to be a consistent way of UV regularization in coordinate space. 
However, in our case a quite sloppy treatment of UV divergences is allowed because we later match our result directly 
to the two-loop expression for the cross section in the threshold region. (See also fej.) 
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the cross section R®Q and the Green function, eq. (Q), leads to another 1/Mq phase space correction 
and using relation <M) the result for the cross section reads 



R QQ = ^N c e 2 Q alm 



H a (a,(3) 



where 



H a (a,f3) = (l 



P 



l + a 



(1 + /? 2 



div 



div 



+ - a 1 Re 
3 



H a {a,f3) 



7 + ln(-»/3) + *(l 



% a ■ 



1 + /3 2 



2(3 



(13) 



(14) 



It should be noted that the term ln(— i (3) in the function H a does not lead to a singular behavior 
of H a in the limit (3 — > because this logarithm is cancelled by a corresponding logarithmic term 
generated by the digamma function in the same limit. It is a remarkable fact that in the framework 



of an expansion in Feynman diagrams the combination 7 + \P(1 



1 a 



2/3 



) is generated entirely by 



diagrams of higher order than the diagrams which produce the explicit term ln(— i/3). (See also [|[.) 
The terms in eq. ( |l~3| ) symbolized by [div] represent divergent and /3-independent contributions. The 
divergences originate from the integration region |x| — ► in TIPT and have to be considered as UV 
divergences which indicate that the electromagnetic current which produces and annihilates the heavy 
quark pair has to be renormalized. This renormalization is usually achieved by the determination 
of the corresponding counterterm via matching to amplitudes calculated in covariant (multi-loop) 
perturbation theory in the framework of QCD ]i~8| , [TjJ. Fortunately this lengthy procedure is not 
necessary in our case because expression (13) can be matched directly to a recent two-loop calculation 



of the fermion-antifermion cross section in the threshold region in the framework of QED [15] which is 
sufficient to determine the C(a 2 ) corrections we are interested in. This "direct matching" procedure || 
is carried out in the formal limit a <C <C 1 for stable quarks, where predictions in the non-relativistic 
effective theory and conventional multi-loop perturbation theory in QCD have to coincide. Expanding 
up to next-to-next-to- leading order in (3 (and including only the 0(a 2 ) contributions with the color 
factors Cp and CfT we are interested in) the two- loop expression (i.e. including Born, one-loop and 
two-loop contributions) for the cross section reads [15] 



K 



QQ 

2 Loop 



N,ej 



+ 



7T 



3vr^ 



7T 



6/? + y/3 2 



1/3 



3C£ + 



5C£tt 2 , 3^ 



24 



+ 



C 2 F In/? )p 



where 



Co 



cl 



1/39 \ 4 , 35 



+ C F T 



4/11 
9 1^2 _1 



(15) 



(16) 



Expanding eq. ( |l3D in the same way and demanding equality to expression (|15|) the divergent contri- 
butions in eq. (|i3|) are unambiguously removed and replaced by constant (i.e. /3-independent) terms. 
The final result for the cross section then reads 



R QQ = ^N c e 2 Q C F a s Im\H a (C F a s ,(3)] ( 1-4C F — + a 2 s C 2 + l C 2 F a 2 Re\ H a {C F a s , (3) 



(17) 



In eq. ( |17D the well known 0(C F a s ) short-distance correction — AC F a s /it [O] is successfully recov- 
ered. It is an interesting fact that the size of the contribution from C2 = —0.24 is an order of 
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Figure 1: The relative 0(a 2 ) corrections to the total non-relativistic tt production cross section 



A 1 (dashed lines), A 2 (dashed-dotted lines) and A 1 + A 2 (solid lines) for a & 
a s = 0.16 (thick lined) as described in the text. 



0.13 (thin lines) and 



magnitude smaller than the one of the function H a . For convenience, C2 will be called C(a 2 ) short- 
distance correction in the following discussion. However, we would like to stress that a unique identifi- 
cation of short-distance and long-distance contributions in the combination a 2 C2 + § C F a 2 Ke[H a ] 
is impossible because such a procedure is cut-off-dependent. In the language of conventional 
perturbation theory (in the number of loops), expression ( |T7| ) resums all contributions of order 
/3(Ci?a s //?) n x [1, Cfci s , (3 2 , Cpa s (3, C F a 2 , Cj?Ta 2 ], n = 0, 1, 2, . . . , 00, in an expansion for small (5. 
Because expression ([17]) is also valid for complex energies it is applicable for ti production, where the 
large top width has to be taken into account^]. 

Although the 0(CaCfO( 2 ) and ©(C^Tr^a 2 ) corrections to the heavy quark-antiquark cross 
section in the threshold region are still unknown, it is instructive to examine the 0(a 2 ) corrections con- 

6 For complex energies there is an ambiguity in the definition of the function H a of order f3 2 which cannot be 
removed by the matching to the two-loop result calculated for stable quarks. For the case of unstable quarks this 
ambiguity amounts to Im[a/3 2 ] ~ CfUsTq / Mq in the total cross section, where Yq is the quark width. For the case of 
ti production (Ft « 1.5 GeV) this ambiguity is of order 0.1% which is beyond the intended accuracy. For bottom and 
charm quarks this ambiguity can be ignored. 
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tained in eq. ( |T7| ) for the case of it production. In Fig. [T] the sum of the relative 0(a 2 ) kinetic energy and 
the phase space corrections, A 1 = (^N c e 2 QCFa s lm[H a ] - RQQ' NK )/RQQ' NIi (dashed lines), the 0{a 2 s ) 
dynamical corrections including the 0{a 2 s ) short-distance contribution, A 2 = a 2 C*2 + § C F a 2 Ke[H a ] 
(dashed-dotted lines), and their sum, A 1 + A 2 (solid lines), are plotted for a s = 0.13 (thin lines) 
and a s = 0.16 (thick lined) in the energy range —10 GeV< E < 10 GeV, where E = y/s — 2Mt and 
Mt = 175 GeV. As mentioned earlier the top decay width is implemented by the naive replacement 
E -> E + i T t , which leads to the relation = [1 - 4M?/(E + iT t + 2M t ) 2 ] 1 / 2 between the velocity 
and the energy E. In Fig. |(a) T t = 1.55 GeV, whereas in Fig. |(b) we have chosen Tt = 0.80 GeV. 
It is striking that the strong energy dependence of A 1 and A 2 around the 15 peak is cancelled in 
their sum leaving a fairly stable correction between 3% and 7% over the whole threshold region^. 
This shows that the IS peak, which is the most important characteristic of the total tt production 
cross section in the threshold region ]7j|, is barely shifted by the 0(a 2 ) corrections determined in this 
work. (See also the near cancellation of the 0(cq) contributions to the energy levels E n in eq. ( |20| ) 
for n = 1.) In particular, the sum A 1 + A 2 is fairly insensitive to variations in the choice of the value 
of indicating that the size of the 0(a 2 ) corrections calculated in this work is not affected by our 
ignorance of a consistent treatment of all finite width effects. The variation of the size of A 1 + A 2 
for different choices of a s further shows that once all 0(a 2 ) corrections are calculated the remaining 
relative theoretical uncertainty for the total cross section can be expected at the level of 1%. Taking 
the size of A 1 + A 2 as an order of magnitude estimate for the sum of all 0(a 2 ) corrections and because 
a consistent 0(a 2 ) analysis has never been accomplished in the framework of potential models for tt 
production we come to the conclusion that an uncertainty of the order 5% percent is contained in all 
predictions for the total tt production cross section based on phenomenological potentials. 
Evaluating formula dlTl) for stable quarks above threshold we obtain 



= -N c e 2 0(l- -0 2 ) 1 V ^J \l-AC F ^ 

r_o,/3>o 2 QH \ 3 y J 1 _ exp (- 2 ^^ZL) \ * vr 



+ a 



Vrel 



C 2 -|c| ( 1 + ln0 + Re*(l-i Cl °~ 



3 ^ V v 2/3 



(18) 



where v re i = 20/(1 + 2 ) is the relativistic relative velocity of the produced quark pair. This verifies 
a suggestion for the form of the relativistic extension of the Sommerfeld expression made in |l5| |. For 
— > this leads to the finite expression 



r>QQ 

^=0,(3=0 



lN c e 2 Q C F a s n[l-4C F ^ + a 2 [c 2 - 2 -C 2 ( In (^) + 7 ) ] } . (19) 

For stable quarks and below threshold formula ([D]) develops narrow resonances at the spin triplet 
(n 3 Si) energy levels (n = 1, 2, . . . , oo)F| 



M Q C 2 F a 2 Mq C F a s 
&n — : — n r 



4 n 2 4 n 4 



11 2 

n 

16 3 



(20) 



7 We would like to mention that A 2 contains contributions of order coming from the /3 2 terms in the function 
H a - These terms represent contributions beyond the intended accuracy and are not included in our analysis. The size 
of these contributions to A 2 does not exceed 0.5% in the considered energy range. 

8 The derivation of the energy levels in eq. ( po| ) from expression (^) is carried out according to eqs. (38)-(40) in m. 
We also would like to note that the O(a^) contributions to the n 3 Si energy levels in eq. ( poj ) are consistent with the 
effective potential V3S1 given in eq. (|ll|). 
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Figure 2: The relative 0{a 2 s ) corrections 5 n for n = 1 (solid line), 2 (dashed line), 3 (dashed-dotted 
line) and oo (dotted line) for values of a s in the range 0.1 < a s < 0.8 as described in the text. 



Parameterizing the resonances in the form (a = 1/137) 
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(21) 



where the M n = 2Mq + E n are the vector resonance masses, we can extract the corrections to the 
leptonic widths (n = 1, 2, . . . , oo), 
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7T 
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(22) 
(23) 



l^n(0)| 2 = MgCpO-l/Sim 3 is the modulus squared of the (unperturbed) Coulomb wave function at 
the origin with the radial quantum number n. It should be noted that the relative C(o; 2 ) correction 
S n remains finite in the limit n — ► oo (called "small binding limit" in [^3] ) , 



drso = lim 5 r , 



C2 - - C F 



In 



+ 7 



(24) 



In Fig. [2] the relative C(a 2 ) correction <5 n is plotted for n = 1 (solid line), 2 (dashed line), 3 (dashed- 
dotted line) and 00 (dotted line) for values of a s in the range 0.1 < a s < 0.8. 

I am grateful to B. Grinstein, J.H. Kiihn, A. Manohar and T. Teubner for useful discussions 
and comments. 
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